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ABSTRACT 

The RPLUS2D code, capable of calculating high 
speed reacting flows, has been adopted to design a com- 
pressible shear layer facility. In order to create react- 
ing shear layers at high convective Mach numbers, hot 
air streams at supersonic speeds, render by converging- 
diverging nozzles, must be provided. A finite rate chem- 
istry model is used to simulate the nozzle flows. Results 
are compared with one-dimensional, chemically equilib- 
rium solutions. Additionally, a two equation turbulence 
model with compressibility effects has been successfully 
incorporated with the RPLUS code. The model has 
been applied to simulate a supersonic shear layer. Pre- 
liminary results show favorable comparisons with the 
experimental data. 


INTRODUCTION 

Recent national interest of trans-atmospheric vehi- 
cle has rekindled the hypersonic research. For this vehi- 
cle, a supersonic combustion ramjet (scramjefc) engine 
was proposed to provide the power. Inside this scramjet 
engine, compressible mixing layers with chemical reac- 
tions are important phenomena. The performance of 
the engine will depend on the supersonic mixing and 
the flame holding of shear layers. 

The behavior of incompressible mixing layers has 
been studied extensively. However, additional study is 
required to understand the compressibility and chemi- 
cal reaction effects at the high speeds. Recognizing the 
importance of the supersonic reacting shear layer, a 
test rig is being designed and built at the NASA Lewis 
center to study the flow physics. This current paper 
reports the design process of the test rig using modern 
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CFD methods and the development of the k — e model 
with compressibility effects for the simulation of high 
speed shear layers. 

The experiment is setup to study the mixing of ox- 
idizer and fuel streams in reacting shear layers for vari- 
ous density, velocity ratios, and Mach number. The pri- 
mary goal is to understand the effects of the compress- 
ibility upon mixing and combustion in a fundamental 
way. Therefore, a two-dimensional shear layer facility is 
designed for its simplicity to quantify the flow physics. 
As shown in Fig. 1, the envisioned facility consists of a 
two-stream wind tunnel with two independent gas sup- 
plies. After passing through flow-management devices 
located upstream, each gas stream expands to its prede- 
termined Mach number by means of a contoured center 
body and tunnel walls. Various combinations of flow 
conditions of high-speed stream and low-speed stream 
allows for the systematic study of mixing and reactions 
of compressible shear layers. 

Currently, there are several test facilities in the 
U.S. to study compressible reacting shear layers. Mun- 
gal et al. 1 designed and built a relatively large facility at 
Stanford University. The test section is 10x5.5x48 cm. 
Vitiated air at supersonic speed is used as high-speed 
(oxidizer) stream while hydrogen mixed with inert gas 
is used as the low-speed (fuel) stream. The facility is a 
blow down type with typical a run time of 30 seconds. 

Dimotakis et al. 2 at Cal. Tech, built a blow- 
down type facility in which premixed volume of hydro- 
gen in nitrogen and fluorine in nitrogen are discharged 
through sonic orifices, maintaining constant mass flow 
rate in each of the free streams. Then both streams ex- 
pand to predetermined Mach numbers and meet at the 
tip of splitter plate. The chemical reaction utilized in 
their experiment is the fastest one known and it reacts 
at room temperature. Therefore, no preheating system 
is needed. Fast chemistry eliminates the issue of the 
finite rate chemistry. However, fluorine is highly cor- 
rosive; tedious preparation and cleanup processes limit 
the availability of the facility. 

Dutton et al. 3 at University of Illinois designed and 


built a blow-down facility. A commercially available 
vitiation furnace fueled by natural gas is used. The 
furnace is operated in a fuel rich condition to elevate 
the total temperature of the hot stream. The combus- 
tion process of methane and air in fuel rich condition 
produces large amounts of hydrogen which is then used 
as fuel for the reacting shear layer. By adopting this 
method, one bypasses the difficult safety procedure for 
handling hydrogen. However, the chemistry is complex 
and can not be easily controlled. 

In our study, the National Aerospace Plane (NASP) 
type application is simulated. In order to obtain de- 
tailed measurement of the whole flow field, the facility 
is envisioned to be a continuous wind tunnel with a 
typical run time of 30 minutes to 1 hour. The shear 
layer consists of a high-speed, supersonic, vitiated air 
(oxidizer) stream merging with a low-speed hydrogen 
(fuel) stream mixed with inert gas (nitrogen) at sonic 
or subsonic speed. For the oxidizer stream, compressed 
air at total pressure about 30 atm is continuously pro- 
vided by a compressor. The total temperature of this 
high-speed stream is elevated to ensure the ignition of 
the shear layer. As illustrated in the Fig. 2, the heat- 
ing process of this stream is composed of two parts, 
namely, a heat exchanger and a vitiation facility. The 
high pressure air first flows through the heat exchanger 
and is heated up to about 860 K, and additional heat 
is needed to further raise its total temperature. A vi- 
tiation heating system in which hydrogen and oxygen 
are burned at stoichiometric ratio introduces hot com- 
bustion products into the air stream and raises its tem- 
perature. In doing so, the chemistry of the vitiation 
process is consistent with that of the shear layer test 
and no additional complexity of chemical kinetics is in- 
troduced into the system. 

Finally, the hot air stream flows through a converg- 
ing diverging nozzle to reach a predetermined super- 
sonic speed and enters into the test section. In the test 
section, the expanded air stream must be about 1000 K 
or above to be able to ignite the hydrogen flow. If Mach 
2.5 is desired in the test section, the total temperature 
of air stream required is about 2000 K. If the heating 
process were to consist only of vitiation, a tremendous 
amount of hydrogen is required, and the oxidizer stream 
is diluted by a large amount of water vapor, which the 
proposed two-stage heating system will help to reduce. 

The flow physics inside the test facility of reacting 
shear layers is complex. It is imperative to know the 
overall thermodynamic properties of various test con- 
ditions to control the experiments. Detailed flow-field 
prediction of each section of the test rig is also desir- 
able. Since there is no chemical reactions in the low 
speed flow stream (fuel stream), the current paper is 
concentrated on the calculations of hot air stream and 


the mixing layers. The objectives are: 

1. Design the nozzle contour to provide a nearly uni- 
form hot air stream. 

2. Predict the flow field of air stream in the nozzle 
with finite rate chemistry. 

3. Compare the results of finite rate chemistry with 
the equilibrium calculations for the hot air stream. 

4. Develop the k — e solver with compressibility effects 
for the simulation of mixing layers. 


THEORETICAL MODEL 


Flow Equations 


The flow equations are taken to be the two- 
dimensional, compressible, Navier-Stokes equations and 
species transport equations. Written in a strong conser- 
vative form, the governing equations can be expressed 
as follows: 

f + |( E - E ') + |(F-F.) = H (1) 

Here x and y are Cartesian coordinates, Q is the depen- 
dent variable, E and F are the convective flux vectors: 
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and H is the source vector: 
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where i = 1, N s — 1. The specific total energy, e, shear 
stress components, r, and heat flux components, q, are 
given as: 
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In the preceding expressions, p is the density, u and 
v are the Cartesian velocity components, p is the pres- 
sure, and e is the specific total energy. The subscript 
i identifies each species, and N s is the total number of 
species. For the ith species, Yj, e t , hi, and u;* are its 
mass fraction, specific internal energy, enthalpy, and 
production rate, respectively. The enthalpy of species i 
is obtained by an integration of C P versus temperature: 


,=/ 
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where C P i is the constant pressure specific heat which is 
expressed as a fourth order polynomial of temperature: 


C pi = C p io + C piX T + C pi2 T 2 + C pi3 T 3 + C pi4 T 4 (10) 


The internal energy of species i can be obtained from h { 
using the ideal gas assumption which is valid for high 
temperature: 

e t - = hi — R{T (11) 


where Ri is the gas constant for species i. The diffusion 
velocity components, u, and are calculated by Fick’s 
law 4 : 
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where 

Am =(1-A',)/^-^- (13) 

t*i ij 

is the effective binary diffusivity of species i in the 
gas mixture, and Xi is the mole fraction of species 
i. Transport properties such as the viscosity and the 
thermal conductivity of each species are considered as 
polynomials of temperature, and those of the mixture 
are calculated based on Wilke’s mixing rule 5 . The bi- 
nary mass diffusivities are calculated using Chapman- 
Enskog theory in conjunction with Lennard-Jones in- 
termolecular potential functions. 6 


Two-Equation Turbulence Model 

In deriving the h and c equations, all the flow prop- 
erties are taken Favre averaged (mass weighted aver- 
aged) except the density (p). The definition of the 
Favre average is <£ = <!>pfp where <f> is any flow property. 
In doing so, the resulting equations are much simpler 
because all the terms associating with the density fluc- 
tuation, e.g. p'u\ were eliminated. Thus, Whenever 
the k — e model is applied, all the variables in the afore- 
mentioned equation are changed to the mass averaged 
variables, and the transport properties are changed to 
the effective, i.e. laminar plus turbulent, properties. 
For convenience, no additional flow equations for tur- 
bulent flows with the two equation model is elaborated 
here. 

The k and e equations in the Cartesian coordinate 
system can be cast into the vector form: 
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where G is the turbulence generation term of the k 
equation and can be expressed as: 


G - fit 


(f + f) + m + °ry 

\dx dy ) \ Oy Ox) 


(16) 


The eddy viscosity fi t required by the flow equations is 
modelled as: 

fit = C^p^- (17) 

The constants used in the k — e model are the standard 
Jones and Launder’s values:' C ^ = 0.09, C\ = 1.44, 
C 2 = 1-92, <Tk = 1.0, and <r t zz 1.3. These constants 
were never altered during the course of this work. 

In the source term of the k equation, the dis- 
sipation term pe is multiplied by a correction factor 
(1 -f aM?) to accommodate the compressibility effect. 
Here M t is the local turbulence Mach number defined 
as M t = \fkja where a is the local speed of sound. The 
constant a in the term is taken as unity. This model 
is developed by Sarkar et al 8 The physical meaning of 
the term is that for high turbulence Mach number (A/ t ) 
flows, the dissipation of the turbulence kinetic energy 
is enhanced by a factor of a A/ 2 . For free shear layers 
at high convective Mach number, this correction factor 
reaches a value of 1.3. Thus, the turbulence intensity 
is greatly reduced due to this compressibility effect. 

In calculating the turbulent free shear layers, the 
inlet boundary conditions for mean velocities and tem- 
perature are specified based on the hyperbolic tangent 
profile with specified initial shear layer thickness. The 
hyperbolic tangent profile is an approximation of the 
self-similar solution for fully developed turbulent free 
shear layers. The inlet transverse velocities are set to 
be zero. The turbulent kinetic energy and dissipation 
are specified according to the local equilibrium assump- 
tion and a algebraic turbulence model : 9 

fU = pll (18) 

where / m = 0. 1 255 and the shear layer thickness, 6 , is 
based on the distance between the two transverse loca- 
tions where U = U\ — 0.1AJ7 and U = U 2 + O.lAt/. 
The dissipation can be related to the local length scale 
which is specified based on the local shear layer thick- 
ness: 

k* 

t = C<— (19) 

where C t = 1.23. Using Eqs. (18) and (19), and the 
equation for the eddy viscosity, i.e, Eq. (17), k and 
e can be readily obtained for the upstream boundary 
conditions. 


Numerical Method 


The numerical solution of Eq. (1) is performed in 
a general, body-fitted coordinate system, (£ , 77 ). Coor- 
dinate transformation of Eq. ( 1 ) gives 

^ + |(E-E.) + |-(F-F.)=H (20, 

where 

Q = hQ 

E = A(&E + £,F) 

F = h(rj x E + rj y F) 

E„ = /»(£ X E V +£j,F„) 
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H = /ill 


in which h is the cell volume. 

The transformed equation, Eq. (20), is solved us- 
ing a time-marching, LU scheme. The LU scheme can 
be obtained by approximately factorizing the left-hand- 
side (LHS) of the equation. Tn time-marching form, the 
implicit upwind difference scheme of Eq. (20) can be 
written as 


[I + Af(D+A~ + D+B~ - D+ 

DJ A + + £>~B + )]AQ = A tRHS 


( 21 ) 


In Eq. (21), At is the time-step. Backward-difference 
operators are denoted by D £ and , and forward- 
difference operators by and D+ . The flux Jaco- 

bians, A + , B + , A~, and B" are constructed such that 
the eigenvalues of matrices are nonnegative and 
those of * — * matrices are nonpositive. The matrix D 
is the source Jacobian. 

The left-hand-side matrix of Eq. (21) is usually too 
large for direct inversion. An approximate-factorization 
procedure is implemented to result in the following LU 
scheme : 
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in which the grid spacing in the general coordinate, A£ 
and A 77 are usually taken to be one. Equation (22) is 
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the generic form for the LU schemes. Its derivation can 
be found in Ref. 10 and will not be repeated here. This 
LU scheme requires inversion of the matrix, 

[i + A<(a + -A- + B + -B--d)] (23. a) 

for the L operator and 

[i + A«(a + - A- + B+-B - )] (23.6) 


for the U operator. Note that the source Jacobian 
appears in the L operator only. This is intention- 
ally designed to minimize computing the source Jaco- 
bian, which is quite expensive for most non-equilibrium 
chemistry models. 

Up to this point, no definition has been made to 
the exact form of the split flux Jacobians. For LU- 
SSOR, the split flux Jacobians are defined as 


A + = 0.5(A + TaI) 
A" = 0.5(A - 7a 1 ) 
B + = 0.5(B + 7 qI) 
B- = 0.5(B - Tnl) 


(24) 


where 7 ^ and 73 are greater than the spectral radii of 
the associated flux Jacobians : 


Ta > "”“<l*Al> (25) 

7b > max( | Aq | ) 

The purpose of constructing split flux- Jacobians by Eq. 
(24) is to make the matrices in Eqs. (23. a) and (23. b) 
diagonal for efficient inversion. Apparently, the eigen- 
values of the split flux- Jacobians are not the character- 
istics speeds of the flow. 

In deriving Eq. (23) for the LU-SSOR the equiv- 
alent 7 is assumed to be frozen. This assumption 
tremendously simplifies the algebra involved, and cur- 
tails the computational operations needed. 

In calculating the RIIS, both the inviscid and vis- 
cous terms are discretized using the central-difference 
approach: 

RHS = D;( E v -E) + D n ( F v - F) + H (26) 

where and D ^ are the central difference operators. 

In solving the k and e equations, the aforemen- 
tioned numerical method, i.e., the LU-SSOR on the 
left hand side and central differencing on the right hand 
side, is used. The solution procedure of the whole equa- 
tion set is decoupled into flow solver and turbulence 
solver. Thus, the turbulence solver stands alone and 
can be easily turn on or off. This arrangement does 


not affect the overall numerical stability due to the fact 
that the feedback from k and c equations to flow equa- 
tions depends on turbulent transport properties only. 
Thus, it is more efficient and convenient to separate 
the solution procedure into two parts. 

The source terms of the k and c equations de- 
mand special treatment. In linearizing the source terms 
for the numerical method, Jacobian matrix is obtained 
through the derivative of the source terms with respect 
to the dependent variables, i.e., pk and pc. Following 
the usual practice, the form of the source terms guar- 
antee a 2 x 2 full matrix for the Jacobian matrix. How- 
ever, special treatment in deriving the Jacobian matrix 
is applied in this work to enhance the numerical stabil- 
ity. In the k equation, c has been replaced by k(c/k) 
where c/k is treated as a constant. A similar method 
is also applied to the source term of c equation. The 
Jacobian matrix obtained is: 

D= (-*(! + ««?) _® jf ) (27) 

Note that off-diagonal terms are eliminated and the di- 
agonal terms are always negative. Thus, the implicit 
part of the source terms of k and e equations behaves 
like a sink which always stabilize the numerical scheme. 


Results and Discussions 


Equilibrium Calculations 

The chemical equilibrium program (CET89 code) 
developed by Gordon and McBride 11 at the NASA 
Lewis Center is adopted for the overall evaluation of air 
stream conditions. By providing initial condition of the 
oxidizer and fuel, the CET89 code performs a free en- 
ergy minimization to obtain chemical composition and 
thermodynamic properties at equilibrium status. Ad- 
ditionally, it also provides isentropic expansion calcula- 
tions for either frozen or equilibrium processes. Since 
we are interested in the thermodynamic properties of 
the air stream after expansion, this special feature of 
the CET89 code, namely, the rocket performance cal- 
culation, is used. 

Table 1 shows a typical result of the equilibrium 
calculation with isentropic expansion. The first part 
of the table shows the initial condition of the flow sys- 
tem: one mole of air at 860 K mixed with hydrogen 
and oxygen at about 300 K. The total pressure is 17 
atm, and hydrogen and oxygen is in stoichiometric ra- 
tio for vitiation heating process. This is a H-P problem 
that specifies the total enthalpy and and total pressure 
of the system and then calculates the adiabatic flame 
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temperature. Iterative calculations adjusting the the 
amount of hydrogen and oxygen added to the system 
yields the desired temperature and Mach number of the 
expanded air flow. As shown in Table 1, for one mole of 
air, 0.22 mole of hydrogen and 0.1 1 mole of oxygen are 
needed for vitiation heating to obtain a hot air stream 
at Mach 2.5 and 1000 K of temperature. 

Equilibrium calculations are used to account for 
the changing species composition during the expansion 
process before the air enters the test section. Inside the 
combustor chamber, minor species such as 1102, NO, 
OH ...etc. exists. During the expansion process, these 
high energy radicals recombine and form more stable 
species such as H2O, O2, and N2; according to calcula- 
tion, no minor species survives. During the expansion, 
specific heat of constant pressure (C P ) and specific heat 
ratio (7) change significantly due to the change of the 
species composition and the diminishing of the vibra- 
tion mode of internal energy. Therefore, the so called 
"real gas effect” must be included in the CFD simula- 
tion of the flow expansion. 

While the equilibrium calculation has provided a 
global understanding of the vitiation and expansion 
processes of the air stream, a more detailed flow field 
solution is invaluable for the understanding of the two 
dimensional effect, boundary layer flow, and the finite 
rate flow chemistry. 


Finite-Rate Calculations 

Before the Navier Stoke calculation, the nozzle 
contour for the flow acceleration is designed. A com- 
mercially available MOC (Method Of Characteristics) 
code based on Rao’s method is adopted in designing 
the nozzle contour. The calculation is based on the 
ideal gas with constant specific heat ratio. The code 
performs wave cancellation calculation to deduce the 
nozzle contour which provides a uniform flow at the 
exit of the nozzle. However, due to the existence of a 
contoured center body in the experiment, the designed 
contour can not be applied to the test rig directly. As 
shown in Fig. 3, the subsonic part of the nozzle must be 
bent. Therefore, CFD code is needed to verify that the 
bent nozzle still provide a uniform flow into the test 
section. Figure 3 shows the geometry of a Mach 2.5 
nozzle for the shear layer facility. The contour of the 
nozzle is designed using the Rao code. The subsonic 
part of the nozzle is bent to accommodate the center 
body which separates the high speed stream and low 
speed stream. 

The rest of the figures shows calculations of the 
RPLUS code. Figure 4 shows the Mach number con- 
tour for the supersonic bent nozzle. The Mach number 


at the exit is about 2.4. The fluid accelerates dramat- 
ically in the throat region and becomes supersonic in 
the diverging part of the nozzle. Due to the high shear 
stresses, the boundary layers in the supersonic part of 
the nozzle are significantly thicker that that in the sub- 
sonic region. Figure 5 shows the Mach number profile 
at the exit of the nozzle (the inlet of the test section of 
shear layers). Except for the existence of boundary lay- 
ers, the flow is quite uniform at the nozzle exit. Thus, 
the bent nozzle did not aggravate the uniform flow de- 
signed by the MOC calculation. The temperature con- 
tour is shown in Fig. 6. Temperature decreases as the 
flow accelerates. Note that the temperature at the noz- 
zle exit is about 1100 K, which compares favorably with 
the equilibrium calculation, theoretically, this temper- 
ature is adequate for the ignition of diluted hydrogen 
flow provided by another stream. Again, strong tem- 
perature gradient occurred inside the boundary layers 
due to high shear stress. The pressure contour is illus- 
trated in Fig. 7 showing that the pressure decreases 
as the flow speeds up. The pressure at the inlet of the 
test section is about 1 atm. This test condition of rel- 
ative lower pressure curtails complex treatment of the 
hydrogen/air combustion chamber. The pressure pre- 
dicted by flow solver also agrees favorably with that of 
the equilibrium calculation. 

Figures 8-11 show the species concentration con- 
tours of H2 , O2, H2O, and OH. The chemical kinetics 
is the dominant factor in determining the species distri- 
bution and kinetics is strongly temperature dependent. 
Thus, the pattern of the species contours resemble that 
of temperature, i.e., very little change in the subsonic 
region and significant recombination in the supersonic 
region. Basically, very minor change occurs for the ma- 
jor species such as 02, H2, and H2O, whereas high en- 
ergy radicals, such as OH, diminished in the flow accel- 
eration process. The mass fraction of oxygen at the exit 
of the nozzle is about 20%, which is very close to the 
composition of air. The vitiation system using oxygen 
and hydrogen does provide desirable oxidizer stream. 

Turbulent Shear Laver Calculations 

The test case of the newly developed k — e solver 
is one of the experimental tests reported by Dutton et 
ai. 3 The test condition is tabulated as Table 2. In this 
case, two streams at different temperatures, densities, 
and Mach numbers merge together to form a nonre- 
acting free shear layer. The convective Mach number 
is 0.45. Figure 12 shows the development of the free 
shear layer. Note that x and y axes are not on 1:1 ra- 
tio for the convenience of observation. The definition 
of the shear layer boundaries is the same that of the 
shear layer thickness. In Fig. 12, the boundaries of 
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the shear layer corresponding to 10 — 90% are drawn. 
Circles are the experimental data by Dutton cl a]. The 
calculated results underpredict the development of the 
shear layer. After the developing region, the boundary 
of the shear layer is almost linear. Incidentally, Fig. 13 
shows the numerical convergence trends of the flow and 
k - e equations. In about 2500 iterations, the residu- 
als drop about 12 orders of magnitude and reach the 
machine accuracy. 

Figure 14 shows the Mach number profile at vari- 
ous axial locations. The velocity gradient in the trans- 
verse direction decreases as the flow goes downstream. 
Figure 15 is the coalesced version of Fig. 14. The 
nondimensionalized y coordinate defined as (y — y c )/6 
is used, where y c is the transverse location at the center 
of the shear layer, and 6 is the shear layer thickness. 
Note that the upstream boundary condition of veloc- 
ities is prescribed according to a hyperbolic tangent 
curve which is an approximation of self-similar solution 
of free shear layers. According to Fig. 15, this self 
similarity of velocity profiles never fail as the flow goes 
downstream. 

Figure 16 shows coalesced turbulence kinetic en- 
ergies at various locations. The turbulence kinetic en- 
ergy is nondimensionalized by A U 2 . This figure clearly 
shows that after the first three stations, i.e., about 150 
mm, the turbulence kinetic energy retains self similar- 
ity. Thus, the developing region for the turbulence is 
about 150 mm. Figure 17 shows the turbulence dissipa- 
tion profiles. As flow goes downstream, the peak values 
of turbulence dissipation at each stations decrease, and 
the turbulence dissipation never reach a fully developed 
condition. However, if the e value is nondimensional- 
ized by A U 3 /S the coalesced profiles appear as shown in 
Fig. 17. Similar characteristic is observed in the eddy 
viscosity profiles (Fig. 18). Figure 19 shows the nondi- 
mensionalized Reynolds stress profiles. Again, the tur- 
bulence Reynolds stress reach the fully developed condi- 
tion at about 150 mm downstream of the splitter plate. 

Figure 20 shows the comparison between the pre- 
dicted fully developed Reynolds stress and the experi- 
mental data reported by Dutton et al. 3 The predicted 
solution underestimated the peak value of the Reynolds 
stress profile by 6 ~ 8%; however, the overall trend 
of the predicted results is correct. Many reasons con- 
tribute to the discrepancy between the predicted result 
and experimental data. Among them, the upstream 
boundary conditions are simplified in the solution pro- 
cedure, i.e., no effort was made to simulate two bound- 
ary layers merging at the tip of the splitter plate. This 
could offset the solution in the developing region and 
shift the fully developed solutions. 


CONCLUDING REMARKS 

CFD application to the design of the facility lies in 
the provision of detailed information of the whole flow 
field. The chemical reacting flow solver (RPLUS2D) 
is applied in the design process of a compressible free 
shear layer facility. However, a one-dimensional equilib- 
rium calculation is an indispensable tool for the overall 
design of the reacting shear system. The thermody- 
namic properties predicted by the equilibrium calcu- 
lation can serve as a yard stick for detailed flow field 
calculations. A converging-diverging nozzle is used to 
provide a supersonic stream as the high speed stream of 
the shear layer facility. Detailed analyses are performed 
to study the flow and chemical characteristics of the 
nozzle flows. The chemistry in the subsonic region of 
the nozzle flow is largely equilibrium and little change 
of chemical composition is observed with all high en- 
ergy radicals recombining in the supersonic region. 

A k — e solver with compressibility effect is added 
into the RPLUS2D code. A nonreacting turbulent 
shear layer is simulated by the newly developed pro- 
gram. Preliminary results show favorable comparison 
with Dutton's experimental data. 
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Table 1: Results of equilibrium calculation for Mach 2.5 ox- 
idizer stream. 


FUEL 

OXIDANT 

OXIDANT 

OXIDANT 

CHEHICAL FORMULA 
H 2.00000 
0 2.00000 
0 2.00000 
N 2.00000 

HOLES 

0.220000 

0.110000 

0.210000 

0.790000 

ENERGY STATE 
KJ/KS-HOl 

0.000 G 

0.000 G 
17872.875 G 
16944.439 8 

TEMP 
DEG K 

290.15 

298.15 
860.00 
860.00 



CHAMBER 

THROAT 

EXIT 

EXIT 

EXIT 

EXIT 

PINF/P 


I. 0000 

1.0145 

15.000 

16.000 

17.000 

18.000 

P, HPA 


3.4451 

1.8986 

0.22967 

0.21532 

0.20265 

0.19139 

T . DEG X 


1933.01 

1707. 

97 

1062.08 

1045.93 

1030.93 

1016.95 

RHO , KG/CU 

H 

5.7642 0 

3.5958 

0 

6.9953-1 

6.6594-1 

6.3588-1 

6.0881-1 

H, KJ/KG 


522.32 

187. 

26 

-704.36 

-725.39 

-744.85 

-762.94 

H. HOI WT 


26.892 

26.895 

26.897 

26.897 

26.897 

26.897 

rp ir i/nrsinn 

1 5239 

1.4596 

1.3037 

1.2997 

1.2959 

1 . 2924 

1 a a mm a 


1.2558 

1.2692 

1.3108 

1.3121 

1 . 3133 

1.3144 

1 UFl .M/SFC 

866.3 

818 

.6 

656.0 

651.3 

646.9 

642.8 

| HACH NUMBER 


0.000 

1.000 

2.388 

2.425 

2.461 

2.494 

HOLE FRACTIONS 








H02 


0.00001 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

HZ 


0.00001 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

*H20 


0.17988 

0.18019 

0.18033 

0.18033 

0.18033 

0.18033 

1 NO 


0.00514 

0.00244 

0.00005 

0.00004 

0.00004 

0.00003 

! N02 


0.00005 

0.00003 

0.00000 

0.00000 

0.00000 

0.00000 

1 


0.64482 

0.64627 

0.64752 

0.64752 

0.64752 

0.64753 

f 0 


0.00003 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

| OH 


' 0.00078 

0.00024 

0.00000 

0.00088 

0.00000 

0.00000 

! 02 


0.16926 

0.17081 

0.17210 

0.17211 

0.17211 

0.17211 


Table 2: Test conditions of the calculated shear layer. 


Item 

Units 

Value 

M r M 2 


1.91, 1.37 

M c 


0.45 

u,. U 2 

[m/s] 

702, 404 

T,, T 2 

[K] 

334,215 

Pi’ Pz 

[Kg/m 3 ] 

0.57, 0.89 

P 

[atm] 

0.55 

8’ 


0.027 

S’/S’. 


0.54 
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Fig. 1: Schematic of free shear layer. 


Fig. 4. Mach number contours of the Mach 2.5 nozzle. 



air, <16 Kg/s 
P* 31 atm 
T-300K 

compressor supply 



Fig. 2: NASA supersonic shear layer test facility. 


Fig. 5: Mach number profile at the nozzle exit. 




Fig. 3: Numerical grid for the Mach 2.5 nozzle (grid size 
90x70). 


Fig. 6: Temperature contours of the Mach 2.5 nozzle. 
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Fig. 7: Pressure contours of the Mach 2.5 nozzle. 



Fig. 8: H 2 mass fraction contours of the Mach 2.5 nozzle. 



0.1937 


Fig. 9: O 2 mass fraction contours of the Mach 2.5 nozzle. 


Fig. 10: H 2 0 mass fraction contours of the Mach 2.5 nozzle 



Fig. li: OH mass fraction contours of the Mach 2.5 nozzle. 



Fig. 15. Coalesced velocity profiles of different axial loca- 
tions. 








Fig. 17: Coalesced turbulence dissipation profiles of differ- 


ent axial locations. Fig. 20: Comparison of Reynolds stresses between experi- 

mental data and predicted results. 
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Fig. 18: Coalesced eddy viscosity profiles of different axial 
locations. 
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